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ABSTRACT 



D , We study the effect of tfie " chromo-electric" dipole moment on tlie electric dipole 

moment (EDM) of the neutron in the two Higgs doublet model. We systematically 
'^ ■ investigate the Weinberg's operator O^g = GGG and the operator Oqg = qaCq, in 

the cases of tan/5 ^ 1, tan/5 <^ 1 and tan/3 ^1. It is shown that Osg gives the 
main contribution to the neutron EDM compared to the other operators, and also 
that the contributions of Oug and O^g cancel out each other. It is pointed out that 
the inclusion of second lightest neutral Higgs scalar adding to the lightest one is of 
essential importance to estimate the neutron EDM. The neutron EDM is considerably 
reduced due to the destructive contribution with each other if the mass difference of 
the two Higgs scalars is of the order O(50GeV). 
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1. Introduction 

The physics of CP violation has attracted one's attention in the circumstance that 
the 5-factory will go on in the near future. In the experiments of the B decay asymme- 
try, the central subject is the test of the standard Kobayashi-Maskawa model(SM)[l] as 
an origin of CP violation. On the other hand, the electric dipole moment (EDM) of the 
neutron is of central importance to probe a new origin of CP violation, because it is 
very small in SM. Begining with the papers of Weinberg [2], there has been considerable 
renewed interest in the neutron EDM induced by CP nonconservation of the neutral 
Higgs sector. Some studies[3,4,5] revealed numerically the importance of the "chromo- 
electric" dipole moment, which arises from the three-gluon operator GGG by found 
Weinberg[2] and the light quark operator qaCq introduced by Gunion and Wyler[3], 
in the neutral Higgs sector. Thus, it is important to study the effect of these operators 
systematically in the model beyond SM. In this paper, we study the contribution of 
above two operators to the neutron EDM in the two Higgs doublets model(THDM)[6]. 
The 3x3 mass matrix of the neutral Higgs scalars is carefully investigated in the 
typical three cases of tan/3 ^ 1, tan/? ~ 1 and tau/S <^ 1, where tan/3 = |f2/fi| with 
Vi = {(t)^i)vac and 01 and 02 couple with down- and up-quark sectors, respectively. In 
these restricted regions of tan/?, the Higgs mass matrix becomes very simple, and then 
we can easily estimate the CP violation parameters of the neutral Higgs sector, which 
lead to the neutron EDM in THDM. We found that the neutron EDM follows mainly 
from the two light neutral Higgs scalar exchanges. Due to the opposite signs of the 
two contributions, the neutron EDM is considerably reduced if the mass difference of 
the two Higgs scalars is in the order of O(50GeV). 

In order to give reliable predictions, one needs the improvement on the accuracy 



of the description of the strong-interaction hadronic effects. Weinberg employed the 
"naive dimensional analyse" (NDA) as developed by Georgi and Manohar[7] in comput- 
ing the effect of the GGG operator on the neutron. However, this method admittedly 
provides at best the order-of-magnitude estimation. Moreover, when gluon fields are 
present, there occurs an indeterminable factor of 47r, which depends on whether one as- 
sociates a factor Qs or -i-Kgs with each gluon field factor in the interaction Lagrangian[3]. 
Recently, Chemtob[8] proposed a systematic approach which gives the hadronic matrix 
elements of the higher-dimension operators involving the gluon fields by using the large 
Nc current-algebra. In his model, the hadronic matrix elements of the operators are 
approximated by the intermediate states with the single nucleon pole and the nucleon 
plus one pion. So, this approach may be a realistic one. We employ his model to get 
the hadronic matrix elements of the relevant operators in this work. The comparision 
between results by this approach and NDA will be discussed briefly in the last section. 

In section 2, the neutral Higgs mass matrix is analyzed and then the magnitudes 
of the CP violation factors ImZj are estimated. In section 3, the formulation of the 
neutron EDM with the hadronic matrix elements of the CP violating operators are 
discussed. Section 4 is devoted to the numerical results of the neutron EDM and some 
remarks and conclusion are given in section 5. 

2. CP violation parameter in THDM 

The simple extension of SM is the one with the two Higgs doublets[6]. This model 
has the possibility of the soft CP violation in the neutral Higgs sector, which does not 
contribute to the flavor changing neutral current in the -B, D and K meson decays. 
Weinberg [9] has given the unitarity bounds for the dimensionless parameters of the CP 



nonconservation in THDM. However, values of these parameters are not always close to 
the Weinberg's bounds[9]. Actually, the CP violation parameter ImZi(this definition is 
given later) is suppressed by l/tan/3 compared with the Weinberg's bound at the large 
tan/9 as pointed out by Barr[5]. Chemtob[10] has predicted CP violation parameters 
by the use of the renormalization approach under the assumption that the coupling 
constants of the Yukawa couplings and self-coupling scalar mesons interactions reach 
infrared fixed points at the electroweak scale. This infrared fixed points approach either 
leads to the large top quark mass rrit ~ 230GeV, which is unfavourable to the recent 
electroweak precision test, or leads to the existence of the unobserved fourth generation 
quarks. Thus, it is difficult to estimate the reliable magnitudes of the CP violation 
parameters lm.Zi{i = 1,2). However, we found that the Higgs mass matrix is simplified 
in the extreme cases of tan/? ^ 1, tan/? ~ 1 and tan/3 ^ 1, in which the CP violation 
parameters are easily calculated. 

The CP violation will occur via the scalar-pseudoscalar interference terms involving 
the imaginary parts of the scalar meson fields normalization constants, Zi, which are 
column vectors in the neutral Higgs scalar vector space, defined in terms of the tree 
level approximation to the two-point function as follows: 

[^(0?0?)., ^(0^0^)., ^(0K)., ^M<pfy 
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where Vi = {(f)^)vac- The CP violation factors \Ta.Z\ are deduced to 
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where u\ denotes the i— th component of the k—th normalized eigenvector of the 
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Higgs mass matrix. Let us estimate u] by studying the symmetric Higgs mass matrix 
M^ whose components are 
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which are derived from the Higgs potential 
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As a phase convension, we take h to be real and 
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Now, the Higgs mass matrix M^ is rotated so as to make the (1,3) (and then (3,1)) 
component zero by the orthogonal matrix Uo as 



COS/9 sin/9 
Uo = I — sin /9 cos (3 \ . 
1 



(6) 



Then, the transformed matrix M'^ = UoM^Uq is given as 

Ml? = 2^icosV + 2^2sin^/3 + 4(?-^)sin2/3cos2/3 , 

M22 = 2{gi + g2 + 2g-2^)sm'^ Pcos"^ P + g' + '^+hcos2(l) , 

M'i = g' + l-hcos2<j), 

M[l = smpcosp[cos2p{2g-2'^ + gi + g2)+gi-g2\ , (7) 

M[l = 0, 

M23 = -h sin 20 , 

in the v"^ = |fip + |f2p unit and the parameter ^ is defined as .^ = ^/|fif2p- This matrix 
cannot be diagonahzed in the analytic form generally, unless special relations among 
the parameters of the mass matrix are satisfied. The parameters are only constrained 
by the positivity condition as follows[10,ll]: 

gi>0 , g2>0 , h<0 , h + g' <0 , g + g' + h> -y/glg^ . (8) 

However, we can simply diagonalize the Higgs mass matrix in the extreme cases of 
tan/9 ^ 1, tan/5 ~ 1 and tan/5 ^ 1. 

At first, we consider the case of tan/3 ^ 1. By retaining the order of cos/3 and by 
setting cos^/5 = 0, sin/5 = 1, the mass matrix becomes 

2(72 2cos;9(e-^-^2) \ 

2 COS I3{^- g- g2) ^' + ^ + /icos20 -/i sin 20 . (9) 

-/i sin 20 g' + '^- hcos2(f)j 

In the limit of cos/3 = 0, this mass matrix is diagonahzed by only rotating on the 

(2-3) plane. However, due to the non-vanishing tiny M^^^ component, this rotation is 

slightly deviated from the (2-3) plane. The orthogonal matrix Ui to diagonalize the 



Higgs mass matrix of Eq. (9) is approximately obtained as: 

/I \ / 1 ecos0 0\ / 1 esin0\ 

Ui ~ COS0 sin0 -ecos0 1 1 ' (10) 

VO -sin0 cos^y V 1/ V-esin0 1 / 

where, neglecting /i cos 20, 

e~^ cosp . 11) 

^ + g'-2g2 

Then, the eigenvectors of M^ in Eq.(l) are 

u^^' = { cos/? — esin/3, — sin/3, 0}, 

■u*-^-* = {sin/3cos0, (cos/5 — esin/3) cos0, — sin0}, (12) 

^(3) _ I sin /3 sin 0, (cos/3 — esin/3) sin0, cos0}, 

with the order of 0(cos^/?) being neglected. The diagonal masses are given as 

M^ = 2g2+0{cos^ (3), MJ = g'+l+h+0{cos^ (3), M| = g'+l-h+Oicos^ f3) , (13) 

in the v'^ unit. The lightest Higgs scalar to yield CP violation is the second Higgs 
scalar with the mass M2 since h is negative and C, is positive. The Higgs scalar with 
Ml does not contribute to CP violation because of mIj = 0. The absolute values of 
g' is expected to be 0(1), but h seems to be small as estimated in some works[ll]. 
Therefore, the masses M2 and M3 may be almost degenerated. Then, CP violation 
is reduced by the cancellation between the two different Higgs exchange contributions 
ImZj and ImZ; since ul u^ and ul u^ (i=l,2) have same magnitudes with opposite 
signs. Thus, it is noted that the lightest single Higgs exchange approximation gives 
miss- leading of CP violation in the case of tan/3 ^ 1. 

(2 3) 
In order to get the magnitudes of ^2 ' , we estimate e, which depends on the value 

of C,- The parameter C, is determined by the charged Higgs mass as follows: 

M±2 = ^^2 _ (^4) 



We have already studied the charged Higgs scalar effect in THDM through the inclusive 
decay B -^ Xs'~f [12], as to which the upper bound of the branching ratio was recently 
given by the CLEO collaboration[13]. We obtained 300GeV for the lower bound of the 
charged Higgs scalar mass in the case of rrit = ISOGeV and m^ = 5GeV. This lower 
bound means .^ > 3. In the limit of the large ^ with retaining other parameters to be 
0(1), e/cos/3 reachs 2 as seen in Eq.(ll). Actually, gi, g2, g and l^f'l are around 1 in 
some numerical studies [11]. Then, if we take M"^ = 400(350)GeV, which corresponds 
to ,^ = 5(4), the value of e/ cos/? becomes 3(4). In the following calculations, we fix to 

(2) 

be e = 3 cos/?. By use of these resulting u] values, we can calculate the CP violation 

(2) 

factors ImZi = ImZl . We show the numerical results together with the Weinberg 
bounds for ImZi and Im.Z2 in Figs. 1(a) and 1(b), where = 7r/4 is taken. Although 
the Weinberg bounds give nothing for these signs, our estimates determine the relative 
sign between IraZi and lm.Z2. For IraZi, our result reachs the Weinberg bound, but 
for Im.Z2 the our calculated value is suppressed compared with the Weinberg bound in 
the order of 1/tan/?. 



Figs. 1(a) ~ 1(f) 



CP violation in the case of tan/? <^ 1 is similar to the one of tan/? ^ 1. By 
retaining the order of sin/? and by setting sin^/? = 0, cos/? = 1, the mass matrix 

becomes 

2j7i _2sin^(^-^-^i) \ 

-2sin/?(^-5f-5fi) g' + ^ + hcos2(j) -/isin20 . (15) 

-/isin20 g' + ^-hcos2(t)) 

Except for replacing g2 with gi and cos /? with — sin /?, the mass matrix is the same one 



as of Eq.(9) in the case of tan/? <^ 1. The eigenvectors are easily obtained as follows: 

ii(i) = {cos/?, -(sin/3 + e'cos/9), 0}, 

■" = { (e' + sin /3) cos 0, cos/3cos0, — sin0}, (16) 

'^ = { (e' + sin /3) cos 0, cos /3 sin 0, cos0}, 

with the order of 0(sin^/3) being neglected. The e' parameter is defines as 

e = — = smp , 17 

^ + g' -2gi + h cos 2(f) ^' ^ ^ 

which is taken to be e' = —3 sin/3 as discussed in Eq.(ll). Taking u\ as the eigenvector 
of the lightest Higgs scalar, we show ImZi and lm.Z2 in Figs. 1(c) and 1(d). For lm.Z2, 
our result reachs the Weinberg bound, while for ImZi the calculated value is suppressed 
from the Weinberg bound in the order of tan f3. The relative sign between ImZi and 
Im.Z2 is just same as in the case of tan/3 ^ 1. 

The last case to consider is of tan /3 ~ 1. Setting cos 2/3 = 0, we get the Higgs mass 

matrix as 

^91 + ^92 + 1-9 1(91-92) 

1(91-92) \gi + \g2 + 9 + 9' + hcos2(j) -/i sin 20 

—/i sin 20 g' + ^ — hcos2' 

The off diagonal components are very small compared to the diagonal ones because 

Qi ~ g2 is suggested by some analyses[ll] and h is small. Then, we get the approximate 

eigenvectors as follows: 

u^^' = { cos/3 — sin/3sin6'i2 cos^23, — sin/3 — cos/3sin6'i2 cos^23, sin ^12 sin 6*23 }, 

u^'^' = { sin /3 cos ^23 + cos /3 sin 6112, cos /3 cos 6*23 — sin /3 sin 6'i2, —sin ^23}, 

M*-^-* = { (sin /3 sin 6*23, cos/3 sin6'23, cos^23}, (19) 



where 

92- 9i 



tan 2^12 



^ — 2g — g' — h cos ! 



4hsm2(j) 2hsm2(f) 2 



tan 2^23 = =^ ^ TT^^ ^^ • (20) 

^i+^2 + 2^-2^ + 4/icos20 Ml -Mi ^ ' 

The Higgs scalar mass Mi is expected to be the heaviest one and the M2 to be the 

hghtest one because of ^ > 3.0 and (?! ~ 5'2 ~ (? ~ l^?'] — 0(1). We estimate the effect 

of CP violation by considering the Higgs scalar with M2 being lightest one and then 

add the effect of the one with M3. Since 612 is expected to be of 0(10^^) [11], we neglect 

terms with sin ^12 in Eq.(19). We can calculate the CP violation parameters ImZj by 

fixing both values of h and M2/M3. We show ImZi and ImZ2 in Figs. 1(e) and 1(f) 

taking M2 = 200GeV, M3 = 250GeV and /i = -0.1. For both ImZ2 and ImZi, the 

calculating values are roughly 1/3 of the Weinberg bounds. The relative sign between 

ImZi and ImZ2 is opposite to the one in the cases of tan/? ^ 1 and tan/? ^ 1. 

3. Formulation of the neutron EDM 

The low energy CP-violating interaction is described by an effective Lagrangian 
LcPi which is generally decomposed into the local composite operators Oi of the quarks 
and gluons fields, 

LcP = Y.C^{M,^^)0,{^Ji) . (21) 

i 

Some authors pointed out [3,8] that the three gluon operator with the dimension six 
and the quark-gluon operator with the dimension five dominate EDM of the neu- 
tron in THDM. So, we study the effect of these two operators on the neutron EDM. 
Various techniques have been developed to estimate the strong-interaction hadronic 
effects[7,8,14]. The simplest one is the NDA approach[7], but it provides at best the 
order-of-magnitude estimates. The systematic technique has been given by Chemtob[8] 

10 



for the case of the operator with the higher- dimension involving the gluon fields. We 
employ his technique to get the hadronic matrix elements of the operators. 

Let us define the following operators: 

0,,{x) = -^qa.^G^'^q , 0,,{x) = -^r^^G^Gl^Gl^ , (22) 

where q denotes u,d or s quark. The QCD corrected coefficients are given by the 
two-loop calculations [2, 3] as follows: 



74 
'23 



C,, = -^^^?^^^,^(/i) [/(.,) tan2/5ImZ,-g(z0cot2/5ImZ,](^^ 



74 
'23 



108 
23 



C3g = ^F7-rM^t)ImZ2 ( ^^ ) , (23) 



where Zt = {mt/niHY a-^id we omitt the upper- indices {k) defined in Eq.(2). The 
function f{zt), g{zt) and h{zt) are the two-loop integral function, which are defined 
in Refs. [4,5,15]. The Csg coefficient is same as Cdg except for the quark mass. In our 
practical calculation, the modification to account for the passage through the h and c 
quarks thresholds involves the replacement 

[gs{M)J ^\g,{M)J \gs{m,)) \gs{m,) J ' ^ ' 

The hadronic matrix elements of the two operators are approximated by the inter- 
mediate states with the single nucleon pole and the nucleon plus one pion. Then, the 
nucleon matrix elements are defined as 

(iV(P)|O,(0)|iV(P)) = Ajl{P)zj,U{P) , 
{N{P')\0,\N{P)n{k)) = BJJ{P')r''U{P) , (25) 

11 



where U{P) is the normahzed nucleoli Dirac spinors with the four momuntum P. Using 
Ai and Bi{i = ug, dg, sg, 3g), the neutron EDM, d'^, are written as 



dl = Y~^ ^ ^'^' + F{g^NN, frin, m^) ^ dBi , (26) 



where yU.„ is the neutron anomalous magnetic moment. The function F{gT^Ni\f,mn,rn.,r) 
was derived by calculating the pion and nucleon loop corrections using the chiral La- 
grangian for the coupled Nn'-f and in given in Appendix A of Ref . [8] . Here, the dimen- 



sional regularization with the standard MS scheme is used for defining the finite parts 
of the divergent integrals. The coefficients Ai and Bi were given by the use of the large 
Nc current algebra and the rjo meson dominance [8]. Then, we have 
A - fn R - 4(m„ + ma)aji 



with ai = -(mso - ms)/(2m, - m„ - m^) ~ -0.28 and F^ = ^2/3Fo = 0.186GeV, 
where fi is defined as 

(77o(g)|0,(0)|0) = /,g^ (28) 

The values of fi were derived by using QCD sum rules as follows[8]: 

/,, = -0.346GeV2 , /g, = -0.842GeV3 , (29) 

where fgg denotes the fiavor singlet coupling. 

Now, we can calculate the neutron EDM. Our inputs parameters are 

Aqcd = 0.26GeV , (m„, m^, m,) = (5.6, 9.9, 200)MeV , fi = m^ , 

M = mt = ISOGeV , g^NN = 13.5 , grj^NN = 0.892 . (30) 

Here, it is useful to comment on the value of fi. As the smaller fi is taken, the QCD 
suppression factor increases, and then, the predicted neutron EDM decreases. Although 
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we do not have the rehable principle to fix /i in the leading-log approximation of QCD, 
we tentatively take n = mn, which leads as{^) = 0.54. If we take n = 0.6GeV, which 
gives rather large as{n) = 0.83, as used by Chemtob[8], our predicted neutron EDM 
will be reduced by a factor 2 ~ 4. 

4. Numerical analyses of the neutron EDM 

We show the numerical results in this section. Since the CP violation parameters 
ImZi have been estimated in the three cases of tan /3, the neutron EDM is predicted for 
each case of tan/?. Since our results are proportional to sin 20, we take the maximal 
case = 7r/4 in showing the numerical results. We show the contribution of the four 
operators Oug, Odg + Ogg and Osg on the neutron EDM, respectively. At first, we show 
the predicted neutron EDM in the region of 5 < tan/? < 10, which corresponds to the 
case of tan/3 ^ 1, in Fig.2(a), where the combined experimental upper bound of the 
neutron EDM[16], 8 x lO^^^e ■ cm, is shown by the horizontal dotted line. The two 
lightest Higgs scalars have been taken into account in our calculations. Defining the two 
lightest Higgs scalar masses to be itlhi and mH2i we fixed tentatively niHi = 200GeV 
and mH2 = 250GeV, which correspond to h = —0.37. In Fig. 2(a), the contributions 
of Oug and O^g are shown multiplying them by the factor 100 because they are very 
small. It is noted that the signs of these two contributions are opposite, and they 
almost cancel each other. The main contribution follows from the one of Odg + Ogg, 
in which the operator Ogg is dominant due to the s-quark mass. This contribution is 
constant versus tan /3 since the tan P dependence of Cdg + Cgg disappears as seen in the 
Eqs.(2) and (23), and then overlapps perfectly to the total EDM(solid line) in Fig. 2(a). 
Thus, the Osg operator dominates the neutron EDM in the case of tan/5 ^ 1. 
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Fig.2(a) 

As the mass difference of tliese two Higgs scalar masses becomes smaller, the neu- 
tron EDM is considerably reduced since the second Higgs scalar exchange contributes 
in the opposite sign to the lightest Higgs scalar one. In Fig. 2(b), we show the the 
predicted neutron EDM versus mHi/mH2 in the case of tan/? = 10 with mui = 200 
and 400GeV. As far as TnHi/iTLH2 > 0.7(|/i| < 0.68), the predicted value lies under the 
experimental upper bound. Thus, it is found that the second lightest Higgs scalar also 
significantly contributes to CP violation. 



Fig.2(b) 



The neutron EDM in the case of tan/? <^ 1 is shown in Fig. 3(a), where we take 
the region of tan/? < 0.25. The contributions of Oug and Osg become very large due to 
the large ImZ2. However, these contribute to the neutron EDM in opposite signs, so 
they almost cancel each other in the region of 1 ^ tan/? > 0.1 as shown in Fig. 3 (a). 
The remaining contribution is the one of Odg + Osg, which is constant versus tan/?. 
In the region of tan/? < 0.1, the cancelation between Oug and O^g is violated and the 
contribution of O^g donimates the neutron EDM in the region of tan/? ^ 0.1. 



Fig.3(a) 



In Fig. 3(b), we show the predicted neutron EDM versus mHi/mH2 in the case of 
tan/? = 0.1. The allowed parameter region of 'mHi/mH2 is obtained by the experiment 
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and is mHi/TnH2 > 0.95(|/i| < 0.07). In othe words, the second lightest Higgs scalar 
mass should be close to the lightest one. We want to note that the predicted EDM 
with rriHi = 400GeV is larger than the one with rriHi = 200GeV in the region of 
f^Hi/f^H2 ^ 0.5. The cancelation between Oug and O^g is violated and the contribution 
of O^g dominates the neutron EDM in the case of irtHi = 400GeV at tan/? = 0.1. 
Thus, one should carefully analyze the signs and magnitudes of the contribution of 
Ougi Odg + Osg and O^g operators in the case of tan/3 <^ 1 since those sensitively 
depend on the values of niHi, itih2 and tan/3. 



Fig.3(b) 



The neutron EDM in the case of tan /3 ~ 1 is shown in Fig. 4 (a). Since the parameter 
h is independent of the Higgs scalar mass difference in contrast to the above two cases, 
we fix /i = —0.1 as a typical value with rrifji = 200GeV and m//2 = 250GeV. The 
contributions of O^g and O^g are shown multiplying them by the factor 10. Similarly to 
be former cases, the signs of these two contributions are opposite and cancel each other, 
and so the dominant contribution is the one of O^g + Osg, which overlapps perfectly to 
the total EDM(sohd line) in Fig.4(a). 



Fig.4(a) 



In Fig. 4(b), the predicted neutron EDM is shown versus mHi/iTLH2 in the case 
of tan/? = 1 with h = —0.05,-0.1. In the region of itlhi/itlh2 = 0.5 ~ 0.9, the 
predicted value is over the experimental upper bound in the case of rriHi = 200GeV 
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with h = —0.1. Thus, the magnitude of \h\ is rigorously restricted by the experimental 
upper bound of the neutron EDM. In both regions of the large and small mHi/iTLH2, 
the predicted neutron EDM is reduced. At iTLHi/mH2 — 1, the cancellation mechanism 
by the second lightest Higgs scalar operates well, while around miii/mH2 — 0, the 
large mass difference of the two Higgs scalars leads to the small 623 as seen in Eq.(20). 



Fig.4(b) 



In all cases of tan (3, the contribution of Odg + Osg dominate the neutron EDM. The 
effects of Oug and O^g seem to become large only in the region of tan (3 <^ 1 although 
these cancel each other considerably. 

5. Conclusion 

We have studied the effects of the four operators O^g, Odg + Ogg and O^g on the 
neutron EDM. The contribution of Ogg dominates over that of other operators except 
for the region of tan/? ^ 1. Moreover, the contributions of Oug and O^g cancel out 
each other due their opposite signs. This qualitative situation does not depend on the 
detail of the strong interaction hadronic model. Actually, in the NDA approximation [7] 
of the hadronic effect, the effects of the two operators almost cancel out although 
the predicted EDM is smaller than ours by a factor 2 ~ 3. Thus, the Weinberg's 
three gluon operator is not a main source of the neutron EDM in THDM. Of course, 
Weinberg's operator may be dominant one in the other models beyond SM, which we 
will investigate elsewhere. The CP violation mainly follows from the two light neutral 
Higgs scalar exchanges. Since these two exchange contributions are of opposite signs, 
the CP violation is considerably reduced if the mass difference of the two Higgs scalars 
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is within the order of O(50GeV). 

Since our results have been shown by taking sin0 = 7r/4, for an arbitrary our 
predicted neutron EDM is simply scaled by the factor sin 20. This factor is expected 
to be of the order one unless is suppressed by an unknown mechanism in THDM. 
Therefore, our results remain unchanged qualitatively. 

Since our predicted neutron EDM lies around the present experimental bound, its 
experimental improvement may reveal the new physics beyond SM. 
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Figure Captions 

Fig.l: The predicted CP violation factors in the case of = 7r/4. The sohd curves 
show (a)ImZi and (b)ImZ2 in tan/? = 5 ~ 10, (c)ImZi and (d)ImZ2 in tan/5 = ~ 0.3, 
(e)ImZi and (f)ImZ2 in tan/5 = 0.8 ~ 1.2 with M2 = 200GeV, M3 = 250GeV and 
h = —0.1. The dashed curves denote the upper bounds given by Weinberg. 

Fig. 2 (a): The predicted neutron EDM in tan/5 = 5 ~ 10 with rriHi = 200GeV and 
^H2 = 250GeV. The dotted curve and dashed curve denote the contribution by Oug 
and O^g, respectively. The contribution of O^g + Ogg overlapps the total neutron EDM 
shown by the solid line. The horizontal dotted line denotes the experimental upper 
bound. 

Fig. 2(b): The mHi/'fnH2 dependence of the neutron EDM in tan/5 = 10 with niHi = 
200,400GeV. 

Fig.3(a): The predicted neutron EDM in tan/5 = ~ 0.3 with mui = 200GeV and 
^_H"2 = 250GeV. The notations are same as in Fig. 2 (a). The dashed horizontal line 
denotes the contribution by Odg + Ogg- 

Fig. 3(b): The 'mHi/iTLH2 dependence of the neutron EDM in tan/5 = 0.1 with rriHi = 
200,400GeV. 

Fig. 4(a): The predicted neutron EDM in tan/5 = 0.8 ~ 1.2 with rriHi = 200GeV and 
fnH2 = 250GeV. The notations are same as in Fig. 2 (a). 

Fig. 4(b): The mHi/mH2 dependence of the neutron EDM in tan/5 = 1 with mni = 
200, 400GeV and h = -0.05, -0.1. 
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